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Abstract. Helioseismic inversions of the Sun's internal 
angular velocity profile show that the rotation changes 
from differential in latitude in the convection zone to al- 
most uniform in the radiative region below. The transi- 
tion occurs in a thin layer, the tachocline, which is the 
seat of strong shear in the vertical direction. In this Note 
we examine whether this rotation profile can lead to shear 
turbulence at the top of the radiation zone. By using the 
standard solar model, we show that such turbulence can 
be generated only in a narrow region 0.695 Rq Js r < 
0.713 Rq at the equator and even in narrower layers at 
higher latitudes. We conclude that the turbulence gener- 
ated by this vertical shear is unlikely to play a significant 
role in the transport of matter and angular momentum, 
and that other mechanisms must be invoked to achieve 
this. 

Key words: Convection, Turbulence, Sun: evolution, 
Sun: interior, Sun: rotation 



1. Introduction 

Helioseismology has revealed that the rotation regime in 
the Sun changes abruptly from latitude dependent in the 
convection zone to almost uniform in the radiation zone 
below. The top of the radiation zone is thus a region of 
strong vertical shear, and for this reason that layer has 
been called the tachocline (Spiegel & Zahn 1992). Acous- 
tic sounding has shown also that some mild macroscopic 
mixing occurs in this region, to smoothcn the composition 
gradient which, with microscopic diffusion only, would be 
somewhat steeper. Moreover, the Li depletion observed in 
solar-type stars can only be explained by some type of 
macroscopic transport between the convection zone and 
the depth where this fragile element is destroyed. 

Among the possible mechanisms which may produce 
such mixing, a plausible one is the shear instability in- 
duced by the vertical rotation profile. In the present Note 
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we shall verify whether this instability is able to play an 
effective role in the observed mixing. 

2. The solar rotation profile 

Through the inversion of LOWL frequency-splitting data 
spanning 2 years, Charbonneau et al. ( 1998| ) have derived 
the following internal solar rotation rate for r > 0.5i?©: 



si(r,e) = n c + 



1 +erf 



rc 



w 



(n s (0)-n o ). (i) 



The observed latidinal differential rotation Q$ is expressed 
as: 



Qs(#) = ^cq + ci cos 2 9 + C2 cos 4 9. 



inHz, f2 eq = 
-67.13 nHz; 



The parameters values are: fic — 2tt x 432. 
2?r x 460.7 nHz, a = -62.69 nHz, c 2 = 
rc = 0.713 Rq is the radius at the bottom of the solar 
convective zone and w — 0.025 Rq the tachocline thick- 
ness. Similar expressions can be found in the litterature 
e.g. Li & Wilson fll99SD , Corbard et al. ( |1998| ), Antia et al. 
( |1998| ). The fact that f2(r, 9) is almost constant at greater 
depth is in conflict with earlier predictions of a faster ro- 
tating core (Pinsonneault et al. 1989; Zahn 1992), and it 
rules out that angular momemtum be transported verti- 
cally through shear-induced turbulence. But this possibil- 
ity still exists in the shear flow of the tachocline. 



3. Conditions for vertical shear instability 

3.1. Richardson criterion 

In a stellar radiation zone, the stable stratification of en- 
tropy tends to inhibit any instability arising from a verti- 
cal shear dV/dz, where V is the amplitude of the horizon- 
tal velocity and z the vertical coordinate. In the adiabatic 
limit, i.e. when radiative losses are negligible, the shear 
instability is suppressed whenever: 



N < 



dV 

dz 



Rir, 



(2) 
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The strength of the stratification is measured by the 
Brunt- Vaisala frequency N: 



N 2 



gs_ 



(V ad - V) , 



(3) 



with the usual notations for the gravity g, the pressure 
scale height Hp, the logarithmic temperature gradients 
V = dhiT/d\nP and S = -(dlnp/d\n.T) P , P, p and 
T being pressure, density and temperature. The critical 
Richardson number Ri c is of order unity for typical flow 
profiles; in the following we shall take Ri c — 1/4. 

However shear turbulence may still arise, and be sus- 
tained, provided that sufficient heat is lost by the turbu- 
lent eddies to lower their buyoancy. Whether this occurs 
is determined by the Peclet number characterizing these 
eddies: 



vl_ 

K' 



in this expression v and I are the velocity and the size of 
turbulent elements and K the thermal diffusivity: 



K 



16a T 1 —IT 
~3 ri npP' 



a is the Stefan constant, k the Rosseland mean opacity 
and Ti the first adiabatic exponent. 

When P c <C 1 the Richardson criterion takes the mod- 



ified form (Dudis 1974; Zahn 1974) 



N 2 < 



dV 

dz 
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(4) 



from which one infers what are the largest turbulent scales 
that can survive in a stratified shear flow. 



3.2. Critical Reynolds number 

However turbulence will be maintained only if the 
turnover rate of the eddies is faster than their viscous 
decay rate: 



v v 



(2 



i/Re c < vl, 



(5) 



v is the viscosity. This critical Reynolds number Re c is 
smaller than the classical one which governs the onset of 
instability in a shear flow whose velocity varies by U be- 
tween two boundaries separated by the distance L, namely 
UL/v. The reason is of course that v < U and I < L. 

For the value of this Re c we turn to an experiment 
performed by Stillinger et al. ( |1983| ), who measured the 
size of the turbulent motions downstream a flow travers- 
ing a grid. The size of the smallest turbulent eddies was 
found to be 15A£k, where I], — (z/ 3 /e) 1 ^ 4 is the classical 
Kolmogorov length, with e being the energy injection rate 
per unit mass. In the inertial cascade the velocities scale 
as v 3 = el, and therefore: 



vl 



,1/3^/3 > £ l/3 (15 .4^)4/3 = (15 . 4) 4/3^ 
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Fig. 1. The shear rate dV/dz with respect to radius be- 
neath the convection zone for the colatitudes 9 = 90° 
(full), 9 = 60° (dashed) and = 30° (dash-dot-dash). 



from which we draw the critical Reynolds number in that 
experiment: Re c sa 40, a value we shall adopt for the 
present purpose. 

The conditions (|) and (§) can be simultaneously ful- 
filled when: 



dV\ 



— >N 2 — 



v Rc r 



dz I 



K Ri c 



(7) 



4. The turbulent region in the Sun 

Let us now examine where the instability condition (0) is 
fulfilled is the Sun. According to (Q), the maximum shear 
rate at the top of the solar radiation zone amounts to 



dV dfl 

— =^rsm9— « 9.10" 

dz dr 



6 s"\ 



as illustrated in Fig. ([[]). 

It is convenient to introduce the non-dimensioned pa- 
rameter 



dV_ 

dz 



K 



according to Eq. (Q), the turbulent shear instability will 
occur if: 



, . , _ Re c 

? <; Qimit = pj- 



160. 



(8) 



Note that £umit, the ratio of Reynolds and Richardson 
critical numbers, is rather uncertain, in particular because 
Ri c depends on the vertical flow profile. 

We have calculated £ as a function of radius and co- 
latitude in the radiative zone beneath the tachocline of a 
standard solar model, in order to delimit the regions where 
the turbulent instability is established. 
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Fig. 2. For the colatitudes 9 = 90° (full), 9 = 60° 
(dashed) and 9 = 30° (dash-dot-dash), stable and unsta- 
ble regions beneath the tachocline located about at the 
limit Rd Cv between the convection zone and the radia- 
tive interior. 

The standard solar model is computed with CESAM 
code (Morel 1997); the evolution includes the pre main- 



sequence phase. Basically, the physics and calibration pa- 
rameters are the same as in Morel et al. ( pOOOQ ; it uses 
OPAL opacities and equation of state, the microscopic 
diffusion coefficients of Michaud & Proffitt ( 1993] ) and the 
recently updated thermonuclear reaction rates of the Eu- 
ropean compilation NACRE (Angulo et al. 1999[ ). The an- 
gular velocity beneath the convection zone is derived from 
Eq. |l|. The dynamical viscosity is taken as in a fully ionized 
hydrogen plasma. 

5. Results and discussion 

Figure || shows log 10 £ as a function of radius for three 
colatitude values. According to the Eq. (||) criterion, the 
turbulent instability may occur at the equator (9 — 90°) 
in the interval 0.695 i?© < r < 0.713 i?©; the turbulent 
layer is thinner at higher latitudes. In other words, the 
unstable region for the vertical shear instability, beneath 
the convection zone, measures barely 2 percent of the solar 
radius and is localized around the equator. 

This instability should not be confused with the shear 
instability arising from the differential rotation in latitude, 
which is not inhibited by the vertical stratification, and 
which may play an important role in smoothing horizontal 
gradients in composition and angular velocity (Spiegel & 
Zahn 1992). But that instability does not contribute to 



the vertical transport. 

We conclude that another physical process is needed to 
transport matter and angular momentum beneath the so- 
lar tachocline, which would operate more efficiently than 
the turbulence generated by the vertical fi-gradient. A 
plausible mechanism is the meridional circulation which 



probably exists in the tachocline (Brun et al. 1999). A 
more powerful process is required also to establish the 
almost uniform angular velocity at greater depth, which 
could be magnetic torquing (Cough & Mclntyre 1998) 
or transport of angular momentum by waves (Schatzman 
1996; Kumar et al. 1999). 
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